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Abstract. Neutrinos may possibly violate the spin-statistics theorem, and hence 
obey Bose statistics or mixed statistics despite having spin half. We find the generalized 
equilibrium distribution function of neutrinos which depends on a single fermi- 
bose parameter, k, and interpolates continuously between the bosonic and fermionic 
distributions when k changes from -I to +1. We consider modification of the Big Bang 
Nucleosynthesis (BBN) in the presence of bosonic or partly bosonic neutrinos. For pure 
bosonic neutrinos the abundances change (in comparison with the usual Fermi-Dirac 
case) by —3.2% for ^He (which is equivalent to a decrease of the effective nnmber of 
neutrinos by AA,, Ri —0.6), -1-2.6% for and —7% for "^Li. These changes provide a 
better fit to the BBN data. Future BBN studies will be able to constrain the fermi-bose 
parameter to k > 0.5, if no deviation from fermionic nature of neutrinos is found. We 
also evaluate the sensitivity of future CMB and LSS observations to the fermi-bose 
parameter. 
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1. Introduction 

Since neutrinos have spin one half they are believed to obey Fermi statistics. A serious 
argument in favor of this belief is an absence of a consistent quantum held theory of 
half-integer spin particles with any other statistics than the Fermi one. For electrons and 
nucleons this issue has been discussed earlier in ref. EElElIll, where it was shown that a 
possible violation of statistics for these particles is strongly restricted by experiment. On 
the other hand, direct experimental checks of Fermi statistics for neutrinos are absent, 
except for a study of the effects of purely bosonic neutrinos on Big Bang Nucleosynthesis 
(BBN) performed 10 years ago jS]. Recently the idea that neutrinos may possess Bose 
or mixed statistics reappeared in ref. jH] where phenomenological analysis of testable 
effects has been presented. A violation of the spin-statistics relation for neutrinos would 
lead to a number of observable effects in cosmology and astrophysics. In particular, 
bosonic neutrinos might compose all or a part of the cold cosmological dark matter 
(through bosonic condensate of neutrinos) and simultaneously provide some hot dark 
matter [H]. A change of neutrino statistics, would have an impact on the evolution of 
supernovae and on the spectra of supernova neutrinos. The presence of a cosmological 
neutrino condensate would enhance contributions of the Z-bursts to the flux of the ultra 
high energy (UHE) cosmic rays and lead to substantial refraction effects for neutrinos 
from remote sources jSj. 

As shown in j2] the change of the neutrino statistics from Fermi-Dirac (F-D) to 
Bose-Einstein (B-E) leads to a decrease of the ^He abundance produced during BBN 
by about 4%. The results of ref. [H] are in qualitative agreement with ref. [S], however, 
according to the estimates of ref. jS] the change of ^He is somewhat weaker. 

Since the double beta decay excludes the possibility of pure bosonic neutrino [Zj, 
but still allows mixed statistics, we will study here an influence of the partially bosonic 
neutrinos on BBN. We introduce the fermi-bose parameter, k, which describes the 
continuous transition from Fermi to Bose distributions in thermal equilibrium as k 
changes from —1 to -|-1, the boundary cases corresponding to purely fermionic and 
bosonic states with any value of k in between allowed. In particular, k = 0 corresponds 
to Boltzmann statistics. We consider possible constraints on the parameter k, from 
present and future BBN and cosmic microwave background (CMB) plus large scale 
structure (LSS) data. 

2. The generalized distribution function 

The form of the kinetic equation for particles with mixed statistics is not immediately 
evident. The statistics dependent factor F[f] under the collision integral e.g. for the 
reaction l-|-2'«-»-3-|-4in the standard case of pure Bose or Fermi statistics has the 
form: 


^ = fl{Pl)f2{P2) (1 ± fsiPs)) (1 ± hiPi)) 

- f3{P3)h{Pi) (1 ± fl{Pl)) (1 ± f2{p2)) , 


( 1 ) 
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where the ± signs correspond to bosons or fermions. 

Derivation of the kinetic equation in general depends on the operator of particle 
number density and the normalization of states with non-zero number of identical 
particles. The quantum creation-annihilation operators of neutrinos obeying mixed 
statistics seems natural to write in the form: 

Ok = a{ cos 7 + sin 7 , 

= afc’^cos7-h afc+sin7, (2) 

where al and are respectively the Fermi and Bose type operators annihilating states 
with momentum k. However, the particle number operator in the standard form, 
Uk = a'^Qk, with Ofc dehned in 0 is not satisfactory if there are several identical particles 
with the same momenta. The emerging problems can be easily observed if one considers 
the matrix element of the particles number operator, Uk, between the multiparticle states 
dehned in the standard way, (a^)'^|0). This state is not an eigenstate of the operator 
Uk- Moreover, it is unclear if there exists an operator of the particle number with the 
appropriate properties, in particular commuting with the free Hamiltonian. It may even 
be that one has to abandon the Hamiltonian approach for description of mixed statistics 
particles. Thus, at this stage we can only make a reasonable guess about the function 
F. 

One simple possibility which has the correct limiting behavior in the case of pure 
statistics is to suggest that the neutrino distribution function in the hnal states of any 
reaction enters F as 

g„ = cos^ 6{1 - f^) + sin^ 6{1 + f^), (3) 

where 6 is the statistics mixing angle and is the neutrino distribution function, which 
in equilibrium interpolates between the F-D and B-E statistics, so that = fpD for 
5 = 0 and = fsE for 5 = ^/2. 

The natural assumption leading to eq. 0 is that the outgoing neutrino contributes 
partly with a Pauli blocking, 1 — /j,, with weight cos^ 5, and partly with Bose 
enhancement, 1 -t- with weight sin^ 5. 

The angle 5 introduced above should somehow be connected to the angle 7 in the 
operator mixture ©• The equality 5 = 7 looks as a reasonable hypothesis. If this is 
true then the angle 7 which parametrizes a violation of neutrino statistics in double 
beta decay can be constrained from BBN and vice versa. If however, the angles 7 and 
5 are different the bounds from BBN and 2/3-decay are not directly related. 


The distribution (jHI) can be rewritten as 


1 

T—1 

III 

(4) 

where 


K = cos^ 5 — sin^ 5 = cos 25 

(5) 


which we call the fermi-bose parameter. Thus, we suggest that the distribution of 
neutrinos in the hnal state enters the factor F in the combination (^. For example, in 
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the case of elastic scattering <->■ 1/2 + I 2 the neutrino distribution functions with 

mixed statistics appear in the collision integral as 

F = - fl{j)2)] [1 - t^fu{k2)] 

- fu{k2)fl{p2)[l - fl{pi)] [1 - t^fv{ki)] . ( 6 ) 

The same factor (1 — nf^) appears (instead of (1 — fu)) in any process involving mixed 
statistics neutrinos. 

Considering the processes of neutrino scattering and production we have found that 
the factor F{f^), and consequently, the collision integral, vanish if 

fu = = bME/T) + k ]-‘ , (7) 

while all other particle distributions are given by the standard equilibrium Bose or Fermi 
functions. That means that given by eq. m is the equilibrium distribution function 
for the case of mixed statistics of neutrinos. 

For K = +1 the function turns into the Fermi distribution, for n = —1 it turns 
into the Bose equilibrium distribution, while for k = 0, i.e. for equal mixture of Fermi 
and Bose statistics, it becomes the Boltzmann one. 

In £g. [U we present the equilibrium distributions dZD for different values of k. 
According to this hgure the distribution becomes softer with an increase of the bosonic 
fraction. The maximum number density shifts to smaller E/T, and the integrated 
number density increases. 

For a positive k (when the fermionic component dominates) the distribution 
function /j, is bounded from above hj < Ijn. 

If K is negative, then for a large lepton asymmetry a neutrino condensation would 
be possible. Indeed, for a given k, the maximum allowed value of the chemical potential 

is 

=m,-Tln(-K), (8) 


this follows from the condition that should be non-negative. In particular, for the 
purely bosonic case we obtain the usual bound p < niy. 

If the neutrino charge asymmetry is so large that could not provide it, the 

neutrinos would form a Bose condensate with the equilibrium distribution function equal 
to 


jieq) — 


+ ^^(k). 


(9) 


exp(£' — + K, 

Here k is the three momentum of the neutrino, and C is a constant whose magnitude 
is determined by the value of the charge asymmetry of the neutrinos. 

An interesting question to address is how unique is the distribution ©. Can one 
introduce other forms of the mixed statistics equilibrium distributions? An alternative 
to (3) could be 


g „ = cos ’< 5(1 - Sfd) + siii ’ i(l + /be) ■ 


( 10 ) 
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Figure 1. The neutrino distribution function for different values of k 
(-1, -0.5, 0,0.5,1). Here x = E/T. 


Such kinetics is equivalent to having two independent neutrino species, bosonic and 
fermionic, and, if equilibrium is established, it would lead to 24/7 neutrino species (see 
next section). To avoid equilibrium the bound on 5 in this case would be quite restrictive. 

3. Effects on big bang nucleosynthesis 

The effect of the change of neutrino statistics on BBN is related to two phenomena. 
First, the energy density of bosons in equilibrium {k = —1) is larger than the energy 
density of fermions by the factor 8/7. If all three neutrinos have B-E statistics, their 
larger energy density would correspond to an increase of the effective number of neutrino 
species by = 3/7. The second, dominant, effect is an increase of the rate of the 
reactions of neutron-to-proton transformations due to the bosonic (anti blocking) factor 
(1 — Because of the larger rate the freezing temperature of these reactions would 
be lower and consequently the frozen u/p-ratio would be smaller. It can be mimicked by 
a decrease of N^. The first effect can be approximately described as AN^ ^ 0.2(1 — n), 
while the second one is noticeably non-linear (exponential). 

We performed the calculations of the abundances of light elements with the Kawano 
nucleosynthesis code jH] which was modified to include the effects of mixed statistics 
of neutrinos described by the distribution 0. This code [Hj is accurate enough for 
calculations of the relative changes of the abundances, while for the absolute values of 
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Figure 2. Upper panel: the change in the effective number of neutrino degrees of 
freedom found from the change of the ^He abundance as a function of the effective 
fermi-bose parameter n. Lower panel: the relative change of the primordial abundances 
of deuterium, helium-4, and lithium-7, as functions of the effective fermi-bose coefficient 
K. We take 77 = 713/717 = 6.5 • 10“^°. 


the abundances we use the results of a more precise modern code jH] • As a reference value 
of the baryon number density we take the WMAP result na. rj = ub/u^ = 6.5 ■ 10 
The results of the computations are shown in hgs. El and El In the upper panel 
of £g. El we present the change of the effective number of neutrino species, dNy, as a 
function of At, which is equivalent to a decrease of the ^He primordial abundance. If the 
neutrinos have a purely bosonic distribution (/t = —1), the effect is similar to having 
^ -0.57. 


However, the effect of modihed statistics cannot be described by a simple change in 
Ny if other light elements are included. In the lower panel of fig. 1 the relative changes 
of the abundances of ^H, "^He, and ^Li with k are shown. As expected the mass fraction 
of ^He drops down and for pure bosonic neutrinos we get the relative decrease about 
^He(At = -1) 


= 1 ) 


- 1 = -3 


( 11 ) 
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This is slightly smaller than the value —3.7% found in [^. 

The amount of goes up with decrease of k. A higher deuterium abundance can 
be explained by a slower conversion rate of deuterium to heavier elements due to fewer 
neutrons and higher expansion rate at BBN epoch when T 0.8 ■ lO^K. In the pure 
bosonic case the increase is about 2.6 %. The ^Li abundance decreases with k, and for 
K = —1 the decrease is about 7%. 



Figure 3. Upper panel: the ratios of abundances of different elements in the cases of 
purely bosonic neutrinos with respect to the standard fermionic case as functions of 
the baryon number density, rj. The vertically hatched (cyan) region shows the WMAP 
2 ct determination of rj. Lower panel: the absolute abundance of ^He as a function of rj 
for the purely bosonic, Boltzmann, and fermionic neutrino distributions, corresponding 
to K = —1, 0, +1 respectively. The two skew hatched regions show the observation of 
primordial helium from ref. m (lower, yellow) and ref. m (upper, magenta), which 
marginally overlap at Icr. 

Let us confront the absolute values of the abundances for partly bosonic neutrinos 
with observational results. We will use the relative changes presented in £g. |21 and 
the central values of the abundances calculated for usual F-D neutrinos in ref. [Oj: 
Yp = 0.2481, X 2 jj/Xh = 2.44 ■ 10“®, and Xth/Xh = 4.9 ■ 10“^°. Other codes 
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may give even slightly larger values for the helium abundance (G. Steigman, private 
communication). We are here (and in the hgures) using the symbol X for the abundance 
of D or Li ratios to H by number, not by mass. At k = —1 we hnd for ^He: Yp = 0.240, 
which makes much better agreement with the value extracted from observations (for a 
review of the latter see e.g. dl)- Different helium observations yield different results, 
e.g., the ref. [12] hnds Y = 0.238 ± 0.002, and ref. [T2j hnds Y = 0.2421 ± 0.0021 (la, 
only statistical error-bars). These results are shown in hgure 2 as the skew hatched 
regions, where the upper (magenta) is the results of ref. |T3| and the lower (yellow) is 
the results of ref. ^2) • Whether the existing helium observations are accurate or slightly 
systematically shifted will be tested with future CMB observations PI. 

The amount of rises at most to X 2 jj/X h = 2.5-10“^, and the agreement between 
BBN and WMAP data remains good, bearing in mind the observational uncertainties. 
Primordial ^Li drops down to Xt^/Xh = 4.55 ■ 10“^*^, again slightly diminishing the 
disagreement between theory and observations. 

We see that at the present time BBN does not exclude even a pure bosonic nature 
of all three neutrinos. Furthermore, the agreement between the value of the baryonic 
mass density, rj, inferred from CMBR and the predicted abundances of '^He, ^H, and ^Li 
becomes even better. In other words, in the standard BBN model there is an indication 
of disagreement between observations of “^He and - they correspond to different values 
of r] with the observed abundances of “^He indicating a smaller value of r] than the one 
given by CMBR, while agrees with CMBR. Motivated by these results the value of 
AA^i, = —0.7 ± 0.35 was suggested in ref. usi. In the case of predominantly bosonic 
neutrinos, as discussed above, the discrepancy between ^H, ^He, and CMBR disappears. 

When the problem of large systematic uncertainties of primordial helium 
determinations will be resolved and the statistical error-bars will dominate the error- 
budget, one can expect to measure N^, with an accuracy at the level of 0.1. This would 
exclude k, < 0.5, if an agreement with the standard BBN values is found. Otherwise, if 
the discrepancy between '^He and remains it may be considered as an indication to 
the mixed statistics of neutrinos. 

Our results change only slightly with variation of the baryon number density t], as 
seen in £g. 2. The upper panel shows the ratio of abundances for purely bosonic to 
purely fermionic neutrinos. The changes are always of the order a few percent for the 
3 abundances considered. The results are in good agreement with ref. [Sj. The verti¬ 
cally hatched (cyan) region shows the 2a WMAP result. The lower panel shows the 
absolute value of the ^He abundance as a function of g, for purely bosonic, Boltzmann, 
and purely fermionic neutrino distribution functions. For other values of k, the result 
will be in between those lines, and can be obtained using curves of £g. |21 The skew 
hatched (yellow) region shows the range of observed values of the helium abundance 
from refs. [nmni, which marginally overlap at la. 

It is well known that CMB can be used to constrain the number of relativistic de¬ 
grees of freedom at the time of photon decoupling (see e.g. jUj). For CMB and LSS, in 
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contrast to BBN, the presence of the bosonic neutrinos increases the number of degrees 
of freedom from 0 to 3/7 = 0.43 as k goes from +1 to —1. The present bounds from 
CMB and LSS data are insensitive to such changes and too weak to constrain the fermi- 
bose parameter (see |TH] for references). The Planck experiment is forecast to constrain 
the relativistic degrees of freedom to the level 5N k, 0.24 HU at la. This means that 
Planck alone will be able to measure k with a precision of about Ak 1. In particular, 
a pure bosonic distribution function for neutrinos can be excluded at about 2a level. 
An “ambitious” future experiment (see details in ref. na. and for earlier predictions 
see pni) will constrain AN to about 0.02, which corresponds to a determination of k 
with precision An 0.1 at la. 

It is known from cosmological considerations that the mass of fermionic hot 
neutrinos are bounded from above by approximately 1 eV, for recent studies see 
e.g. inii23. This bound is applicable to any hot dark matter particle, independently of 
their statistics, which have the same number density as that of neutrinos. For particles 
which have different number densities, or freeze out at different times, this number 
changes somewhat IZQ. Since large scale structure basically constrains the quantity 

nX = ^ ^ 

93eV nth 

where nth is the number density of a thermal fermionic neutrino, then bosonic neutrinos 
will have their masses constrained a factor 4/3 weaker than fermionic neutrinos. We 
thus see that cosmological probes of neutrino masses remain roughly as strong as always, 
in comparison to terrestrial tritium or double beta decay experiments. 

In addition to the thermal neutrino component bosonic neutrinos might condense in 
the early universe and have much larger number density than thermal relics, as argued in 
ref. jH] . However, the cosmological upper bound on their mass would remain practically 
the same because the latter is valid only for hot thermal relics which suppress structure 
formation at small scales. 

The higher number density of bosonic neutrinos will imply a marginally later freeze 
out, and hence a larger sharing of the entropy from the annihilating electrons than in 
the standard scenario |2Sl El ESI- This effect is, however, very small. More general 
non-thermal neutrino spectra and the effect on freeze out is studied in ng. 

Let us comment on a possibility of large neutrino condensate in the case of partly 
bosonic neutrinos. Such a condensate would contribute as cold dark matter in the 
Universe (and show up in the CMB and LSS analysis). The condensate is formed 
when the lepton asymmetry is larger than that which could be ensured by the maximal 
possible chemical potential. In the case of pure bosonic neutrinos the chemical potential 
is restricted by the neutrino mass and is therefore negligible, especially at the BBN 
epoch. In contrast, in the case of partly bosonic neutrinos the maximal potential given 
by (jHI), or ^ = n/T = rUy/T — ln(— k), can be large. So, in the case of partly bosonic 
neutrinos, the formation of the condensate would imply a large chemical potential, 
which could destroy the excellent agreement with BBN. Due to mixing between the 
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active neutrinos the chemical potentials should be equal for all 3 neutrino species at 
the time of BBN pHl 177] . Then using the strong bound on the leptonic asymmetry 
in the electron neutrinos we hnd the bound k < —0.9: For k > —0.9 such a chemical 
potential for the electron neutrinos will signihcantly underproduce helium, leading to a 
disagreement with observations. For negative chemical potential ^He would be strongly 
overproduced leading to essentially the same bound, n < —0.9. 

That is, neutrinos should be almost purely bosonic to produce the condensate 
and satisfy the BBN bound. On the other hand, almost purely bosonic neutrinos are 
excluded (disfavored) by the double beta decay j7]: the mixing angle 6 at the level 
sin6' ~ 0.8 is still allowed, however, the angle 9 is not necessarily equal to 5 introduced 
above. Notice further that the relation of k with the fermi-bose parameter relevant for 
the fdjd -decay is not clear, as discussed in section 2. 

Anyway, an improvement of the BBN bound on k can exclude a possibility of the 
neutrino condensate which might contribute substantially to the cold dark matter in 
the Universe. 

4. Conclusions 

We hnd the equilibrium distribution function for partially bosonic neutrinos which 
depends on a single fermi-bose parameter, k. The change of this parameter from -|-1 to 
— 1 corresponds to a continuous transition between Fermi and Bose distributions. 

We have considered the inhuence of bosonic or partially bosonic neutrinos on BBN. 
In the extreme case of completely or predominantly bosonic neutrinos the primordial 
abundances change in comparison with the usual F-D cases in the following way: '^He 
decreases by 3.2%, increases by 2.6% and ’^Li decreases by 7%. The agreement 
between theory and observations becomes noticeably better. 

Future determinations of ^He will allow to exclude values of fermi-bose parameter 
K < 0.5, if agreement with the standard case is found. The BBN bounds on k, can be 
compatible with those obtained from the analysis of two neutrino double beta decay [7j. 

Future CMB-I-LSS observations can constrain or observe this parameter, possibly to 
the level Ak ~ 0.1, potentially providing indications of a violation of the Pauli exclusion 
principle. 
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